Astronomy &: Astrophysics manuscript no. nis0562 
(DOI: will be inserted by hand later) 



February 2, 2008 



Data reduction methods for single-mode optical interferometry 

Application to the VLTI two- telescopes beam combiner VINCI 

p. Kervella^''^, D. Segransan^ and V. Coude du Foresto^ 



o 
o 

(N 

1—5 

oo 
(N 



> 

in 

(N 

\o 
o 

o 

Oh 

I 

o 



% 



^ LESIA, UMR 8109, Observatoire de Paris-Meudon, 5 pi. Jules Janssen, F-92195 Meudon Cedex, France 

e-mail: Pierre.Kervella@obspm.fr 
■^ Observatoire de Geneve, 51 ch. des Maillettes, CH-1290 Sauverny, Switzerland 
^ European Southern Observatory, Alonso de Cordova 3107, Casilla 19001, Santiago 19, Chile 

Received 23 October 2003 ; Accepted 16 June 2004 

Abstract. The interferometric data processing methods that we describe in this paper use a number of innovative 
techniques. In particular, the implementation of the wavelet transform allows us to obtain a good immunity of the 
fringe processing to false detections and large amplitude perturbations by the atmospheric piston effect, through 
a careful, automated selection of the interferograms. To demonstrate the data reduction procedure, we describe 
the processing and calibration of a sample of stellar data from the VINCI beam combiner. Starting from the raw 
data, we derive the angular diameter of the dwarf star aCeuA. Although these methods have been developed 
specifically for VINCI, they are easily applicable to other single-mode beam combiners, and to spectrally dispersed 
fringes. 
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1. Introduction 

Although interferometric techniques are now used rou- 
tinely around the world, the processing of interferometric 
data is still the subject of active research. In particular, 
the corruption of the interferometric fringes by the turbu- 
lent atmosphere is currently the most critical limitation 
to the precision of the ground-based interferometric mea- 
surements. 

Installed at the Very Large Telescope Interferometer 
(VLTI), the VINCI instrument coherently combines the 
infrared light coming from two telescopes in the infrared 
H and K bands. The first fringes were obtained in March 
2001 with the VLTI Test Siderostats, and in October 2001 
with the 8m Unit Telescopes (UTs). To reduce the large 
quantity of data produced by this instrument, we have de- 
veloped innovative interferometric data analysis methods, 
using in particular the wavelet transform. We have app- 
plied them successfully to a broad range of interferometric 
observations obtained with very different configurations of 
the VLTI (0.35 m siderostats, 8 m Unit Telescopes, 16 m 
to 140m baselines, K and H band observations). 

Since the first fringes of VINCI, more than 800 nights 
of observations have been conducted with this instrument. 
This allowed us to record a large number of individual 
star observations, under extremely different atmospheric 



and instrumental conditions. The data processing meth- 
ods that are described in the present paper were success- 
fully applied to all these configurations, and consistently 
provided reliable and precise results. This gives us good 
confidence that they are efficient and robust, and can be 
generalized to other interferometric instruments. 
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Our goal in this paper is to give a step by step de- 
scription of the processing of the VINCI data, from the 
raw data to the calibrated visibility. To illustrate this pro- 
cessing, we selected from the commissioning data a series 
of observations of a bright star and its calibrator, a Cen A 
and 6'Cen respectively (Sect.O). A complete overview of 
the data analysis work flow is presented in Fig. ^ It can 
be used as a reference to follow the logical progression of 
this paper. The photometric calibration of the interfero- 
grams is described in Sect. 0] and [S] The criteria for the 
selection of the interferograms are detailed in Sect. El and 
the computation of the visibility values and associated er- 
rors is given in Sect. [3 A number of quality controls is 
applied to the reduced data; they are described in Sect.|S| 
The calibration of the visibility is illustrated in Sect.O We 
demonstrate in particular the computation of the statisti- 
cal and systematic errors on the visibility values. 
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Fig. 1. Overview of the VINCI data analysis work flow. 
The shaded area delimits the processing executed auto- 
matically by the instrument data pipeline. The hatched 
area (lower right) covers the astrophysical interpretation 
of the measured visibility, not adressed in the present pa- 
per. 



2. Instrument description 

2.1. The VLT Interferometer and VINCI 

The Very Large Telescope Interferometer (VLTI, 
Glindemann et al. I2U(J0; Glindemann et al. l2Ul)3al 
Glindemann et al. 12003 bl Scholler et al. IM75|l has been 
operated by the European Southern Observatory on top 
of the Cerro Paranal, in Northern Chile since March 2001. 
In its current state of completion, the light coming from 
two telescopes can be combined coherently in VINCI, 
the VLT Interferometer Commissioning Instrument 



Fig. 2. View of the VINCI instrument installed in the 
VLTI interferometric laboratory. The MONA beam com- 
biner is the visible above the center of the image (white 
box), with its optical fiber inputs and outputs. The beams 
coming from the VLTI Delay Lines enter the optical table 
from the bottom of the picture. Figure not available on 
the astro-ph version. 

Fig. 3. Principle of the VINCI instrument. Figure not 
available on the astro-ph version. 



(Fig. El, or in the MIDI instrument (Leinert et al. 
IMTUj) . In December 2002, MIDI obtained its first fringes 
at A = 8.7 /im between the two 8m Unit Telescopes 
Antu (UTl) and Melipal (UT3). Another instrument, 
AMBER (Petrov et al. Onnn| will soon allow the simul- 
taneous recombination of three telescope beams (its first 
observations are scheduled for 2004). 

A detailed description of the VINCI instrument, in- 
cluding its hardware and software design, can be found in 
Kervella et al. IPUHJI . Fig.OI shows the setup of VINCI. 
The two beams enter the instrument from the bottom 
of the figure, after having been delayed by two optical 
delay lines fDerie l?()0()(l . Once the stellar light from the 
two telescopes has been injected into the optical fibers 
(injections A and B), it is recombined in the MONA 
triple coupler. VINCI is based on the same principle as 
the FLUOR instrument (Coudc du Forcsto et al. I1998|l . 
and recombines the light through single-mode fluoride 
glass optical fibers that are optimized for K band op- 
eration (A = 2.0 — 2.4 /im). It uses in general a regu- 
lar K band filter, but can also observe in the H band 
(A = 1.4 — 1.8 /im) using an integrated optics beam com- 
biner (Berger et al. l2()()l|l . The first observations with this 
new generation coupler installed at the VLTI focus have 
given promising results (Kervella et al. I2003al Kern et 

aL onnkji . 

2.2. Beam combination 

The central element of VINCI is its optical correlator 
(MONA), based on single-mode fluoride glass fibers and 
couplers. It was designed and built specifically for VINCI 
by the company Le Verre Fluore (France). The waveguides 
are used to filter out the spatial modes of the atmospheric 
turbulence. In the couplers, the fiber cores are brought 
very close to each other (a few /im) and the two electric 
fields interfere directly with each other by evanescent cou- 
pling of the electromagnetic waves. Motorized polarization 
controllers allow the matching of the beam polarizations, 
in order to obtain the best possible interferometric trans- 
fer function. 

The general principle of the MONA box is shown in 
Fig.^] MONA contains three couplers: two side couplers 
(that provide two photometric outputs Pa and Pb to mon- 
itor the efficiency of the stellar light injection in the opti- 
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Fig. 4. Principle of the MONA beam combiner. Figure 
not available on the astro-ph version. 



cal fibers) and a central coupler that is used for the beam 
combination. The latter provides two complementary in- 
terferometric outputs /i and 12- The four output fibers 
are eventually arranged on a 125 //m square and imaged 
onto an infrared camera (LISA), built around a HAWAII 
detector. Only four small windows, of one pixel each, are 
read from the detector to increase the frame frequency 
and reduce the readout noise. 

The Optical Path Difference (OPD) between the two 
beams is modulated by a mirror mounted on a piezo trans- 
lator. This modulation allows one to sweep through the 
interference fringes (at zero OPD), that appear as tempo- 
ral modulations of the /i and I2 intensities on the detec- 
tor. While the OPD is scanned, the four output signals 
are sampled at a few kHz. The four resulting time se- 
quence signals (two photometric and two interferometric) 
are then available for processing. The interferogram acqui- 
sition rate can be set between 0.1 Hz (faint objects) and 
20 Hz (bright targets). 

2.3. Coherencing 

During the observations, a simple fringe packet centroid- 
ing algorithm is applied in near real-time to the raw data. 
The fringe packet center is localized with a precision of 
about one fringe (2 fim) after each scan and the resulting 
error is fed back to the VLTI delay lines as an OPD offset. 
This capability, called fringe coherencing, ensures that the 
residual OPD is less than a coherence length despite possi- 
ble instrumental drifts. Still, the correction rate (once per 
scan, i.e. a few Hz) is too slow to remove the differential 
piston mode of the turbulence. A fringe tracking unit is 
anticipated for the VLTI (FINITO, Gai et al. 2003) that 
will remove the differential piston and stabilize the inter- 
ference pattern at the sub- fringe level (fringe cophasing), 
thus enabling longer integration times for the scientific 
instrument. 



3. The selected sample data sets 

3.1. Targets 

To illustrate the processing of the VINCI data on rep- 
resentative files, we have chosen two series of interfero- 
grams obtained respectively on a calibrator star, 6 Cen, 
and a target of scientific interest, a Cen A, on the interme- 
diate length EO-Gl baseline (66m ground length). 6'Cen 
was chosen from the Cohen et al. I|1909|l catalogue. These 
authors compiled a grid of calibrator stars whose angu- 
lar diameter is typically known with a relative precision 
better than 1%. Borde et al. H2002(l recently revised this 
catalogue specifically for its application to long baseline 
interferometry. 



Table 1. Relevant parameters of 9 Cen and a Cen A. 





eCen 


a Cen A 




HD 123139 


HD 128620 


mv 


2.1 


-0.0 


niK 


-0.1 


-1.5 


Spectral Type 


KOIIIb 


G2V 


Teff (K)'' 


4980 


5790 


logs" 


2.75 


4.32 


[Fe/H]^ 


0.03 


0.20 


6\ju (mas)'' 


5.305 ± 0.020 


8.314 ± 0.016 



Teff, logg and [Fe/H] from Cayrel de Strobel et al. (l^nUTl 
and Morel et al. 12000L respectively for Cen and a Cen A. 
Measured diameters from Kervella et al. Il2003bt . 



The observations of a Cen A and 9 Cen discussed here 
were carried out with the two 0.35 m test siderostats of 
the VLTI. Both stars are bright, but 9 Cen is significantly 
smaller than a Cen A, therefore its visibility is higher. The 
relevant properties of 9 Cen and a Cen are reported in 
Table ^ The angular diameter of a Cen A was measured 
for the first time by Kervella et al. H2003b|l . based on a 
series of observations that include the two data sets dis- 
cussed here. 

The file names and characteristics of the two se- 
lected data sets are given in Table|21 to allow the in- 
terested reader to retrieve them from the ESO Archive 
( http: //archive .eso. org/ ) . 

3.2. Acquisition parameters and data structure 

Following the standard procedure used with VINCI, a se- 
ries of 500 interferograms was obtained on each object. 
The two data sets were taken on July 15, 2002, starting 
at UT times 01:32:32 for 9 Cen, and 02:33:09 for a Cen A. 
The piezo mirror scanning speed was set to 650 fiia/s, giv- 
ing a fringe frequency of 297 Hz. This intermediate speed 
is used commonly for the operation of VINCI with the 
VLTI Test Siderostats. The LISA camera frequency was 
set to 1.5 kHz in order to obtain a sampling of 5 points 
per fringe. The choice of the scanning speed (hence the 
sampling rate of the camera) is the result of a compro- 
mise between the photometric SNR and the phase per- 
turbations of the atmosphere (dominant at low scanning 
speed) . The VINCI instrument allows one to scan up to a 
fringe frequency of 680 Hz (camera frequency of 3.4 kHz). 
This extreme speed is useful in the case of observations 
with the 8 m Unit Telescopes to reduce the influence of 
the photometric fluctuations on the interference fringes 
(multi-speckle regime). Fig.Oshows the raw signals of one 
interferogram obtained on 9 Cen. This is the second scan 
in the series of 500, and it is of average quality in terms 
of injected flux stability. The photometric fluctuations are 
clearly visible in all four channels, while the interference 
fringes are located close to the center of the scan. The 
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Optical Path Length (microns) 

Fig. 5. The raw signals /i, /2, Pa and Pb, for one interferogram obtained on 9Cen. The original signals have been 
translated vertically respectively by +90, +40, +40 and ADUs for clarity. 



fringes are naturally in phase opposition between the two 
channels /i and 12- 

Each star observation consists of four files (batches), 
that each contain a series of acquisitions (scans) of the four 
signals coming out of MONA, with four different configu- 
rations of the instrument. The first three batches are used 
for the calibration of the camera background and instru- 
ment transmission, and usually contain 100 scans. The 
fourth batch contains the fringes. They are recorded in 
the following sequence: 

1. Off-source: the two injection parabolas of VINCI are 
displaced in order to remove the star image from the 
single-mode fiber head. 

2. Beam A: the A injection parabola is brought back to 
the position where the star light is injected in the op- 
tical fiber, while B is still off. 

3. Beam B: symmetrically, only the beam B is injected in 
the MONA beam combiner, while A is off. The Beam 
A and Beam, B sequences are used in the computation 
of the K matrix (Sect.^+J- 

4. On-source: both beams are injected into MONA, and 
interference can occur. This last series usually contains 
500 scans and is used to compute the squared coher- 
ence factor ^^ (the interferometric observable, defined 
in Sect.ES- 



4. Preliminary steps 

4.1. Computation of the k matrix 

To properly calibrate the photometric fluctuations of the 
interferometric signals Ii and I2 using the two photomet- 
ric outputs P^ and P^, it is necessary to know precisely 
the coefficients linking the intensities of these four out- 
puts. The relationships between the four channels can be 
approximated, within a very good precision (Coude du 



Table 2. Sample data sets. N is the number of scans. 



File name 



Target Type N 



VINCI. 2002-07-15T01:30:12.042 
VINCI. 2002-07-15T01:30:53. 273 
VINCI. 2002-07-15T01:31:36. 505 
VINCI. 2002-07-15T01:32:31. 645 

VINCI. 2002-07-15T02:30:49. 318 
VINCI. 2002-07-15T02:31:30.046 
VINCI. 2002-07-15T02:32:10. 840 
VINCI. 2002-07-15T02:33:08. 661 



eCen 


Off 


100 


eCen 


A 


100 


6'Cen 


B 


100 


6* Can 


On 


500 


aCenA 


Off 


100 


aCenA 


A 


100 


aCenA 


B 


100 


aCenA 


On 


500 



Foresto et al. I1997|l , by the following expressions: 



h = Ki^A Pa + ni.B Pb 
h = K2,A Pa + K2,B Pb 



(1) 
(2) 



The K coefficients correspond to the differential gains be- 
tween the four channels of VINCI. They include the cou- 
pling ratios of the MONA box, the coupling efficiency of 
each fiber to the physical pixels of the infrared camera, 
and the differential quantum efficiency between these pix- 
els. Due to the chromatic transmission of the couplers, 
the color of the observed object plays a role in the k co- 
efficients, and they also tend to evolve due to the slow 
motion of the fiber spots on the LISA detector. It is there- 
fore necessary to measure these coefficients (the k matrix) 
immediately before each star observation. Each pair of 
K coefficients is computed simultaneously using a classi- 
cal x^ minimization algorithm with two variable parame- 
ters. The errors on the estimation of the k coefficients are 
derived from the residual dispersion of the measurement 
points around the linear model. 

Table |21 gives the k coefficients derived for the 6 Cen 
and a Cen observations. The small differences between the 
K values for the two stars may come from the slightly dif- 
ferent colors of these objects, or from a small variation of 
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Table 3. k coefBcients measured for 0Cen and aCenA. 
6 Cen a Cen A 



200 



KiA 0.7569 ± 0.0061 

KiB 4.1231 ±0.0160 

«:2A 1.3089 ± 0.0054 

«:2S 2.4855 ± 0.0143 



0.7576 ± 0.0037 
4.1877 ±0.0120 
1.2790 ± 0.0044 
2.4735 ± 0.0061 



the alignment of the output spots on the LISA infrared 
camera pixels between the two observations (they are sep- 
arated in time by Ai ~ 1 h) . 

Ideally, the k coefhcients should be balanced between 
the four outputs in order to maximize the efficiency of the 
interference, and simultaneously to give high SNR pho- 
tometric signals for the calibration of the interferograms. 
The observed inbalance (that can reach up to a factor 5 
in the selected sample batches) is due to the fact that 
the unexpectedly fast aging of the three optical couplers 
in the MONA box has increased significantly their sensi- 
tivity to temperature. This effect cannot be corrected on 
the coupler itself, and causes a slow (timescale of months) 
but large amplitude evolution of the k matrix. Due to the 
very different time scales of these variations (months) and 
of the scientific observations (hours), this sensitivity is ex- 
pected to have no significant impact on the observations 
other than a uniform and moderate reduction of the qual- 
ity of the LISA signals. 



4.2. Fringe localization 

The first step of the processing is to trim the long interfer- 
ogram to restrict it to a shorter segment, where the fringes 
are centered. The detection of the fringes is then achieved 
with the Quicklook signal QL, that is computed using the 
simple formula: 



QL ^ Ii - a I2 



where a is given by 



^N 



Y:Uh{^)h{^ 



^N 



T:^.iii^ 



(3) 



(4) 



where N is the number of samples of the raw signals /i and 
I2 ■ This operation attenuates the photometric fluctuations 
and increases the SNR of the fringes. Once the fringes 
are localized precisely by detecting the maximum of the 
wavelet power spectral density (WPSD) of the QL, the 
four raw signals are trimmed accordingly. If the fringes 
have been found too close to the edge of the interferogram, 
the scan is discarded to avoid any bias of the data. In 
addition, if a large amplitude jump of the position of the 
fringe packet is detected between two consecutive scans 
(more than 20 /im), a strong piston effect is suspected and 
the scan is rejected (see also Sect.EJ- 




40 60 80 

Optical path length (microns) 

Fig. 6. Photometric calibration of the /i signal. The raw 
/i signal is the black line in the upper part of the plot, 
the photometric calibration signal kiaPa ~ hibPb is the 
superimposed grey line, and the result of the subtraction 
is the lower curve. 



5. Photometric calibration 

5.1. General principle 

The photometric calibration of the interferograms pro- 
duced by VINCI is achieved using the two photometric 
control signals Pa and Ps and the K-matrix. The cali- 
bration is computed separately for the Ii and I2 chan- 
nels using the following formulae (see Coude du Foresto 
et al. l1997 for their derivation) : 



1 



'lea/ 



h - i^i.aPa - i^i,bPe 



2 ^Ki^A Ki^B WPa -Psjwicncr 

1 I2- K2aPa - H2,bPb 



2 ^H2,A K2,B WPa Pb] 



Wiener 



(5) 
(6) 



The subscript "Wiener" designates optimally filtered sig- 
nals fSect. lST^ . This process allows one first to subtract 
the photometric fluctuations that were introduced in the 
interferometric channels by the turbulent atmosphere (cal- 
ibration), and then to normalize the resulting signals to 
the geometrical mean of the two photometric channels 
P = VPa Pb- As an example of calibration, the subtrac- 
tion of kiaPa + hi,bPb from the original /i is presented 
in Fig.lHl The Wiener filtering of P, essential to avoid nu- 
merical instabilities, is described in the next paragraph 
(Sect. 1^21) ■ After the normalization, /i cai and I2 cai are 
apodized at their extremities, to prevent any edge effect 
during the numerical wavelet transform. 

5.2. Low pass filtering of photometric signals 

The normalization by the P signal is a critical step of the 
calibration. If P presents too low values ("zero crossing"), 
the divisions of Eq.O and El will amplify the noise of the 
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Fig. 7. Average power spectral density of the Pa and Pb 

— 2 
signals as measured on 9 Cen. From top to botton: Pa , 

2 N '~2 

noise model Na (upper dotted line), Pb , noise model 

2 

Nb (lower dotted line). 



numerator. This is the reason why the Pa and Pb signals 
have to be filtered, to improve their SNRs. This is achieved 
using Wiener filters, that allow one to optimally filter the 
raw signal and to reject the detector noise. They are com- 
puted from the average power spectral density (PSD) of 
the photometric fluctuations of Pa and Pb using only the 
on-source spectra. We use the classical definition of the 
Wiener filter Wx as computed from the signal Px and the 
noise Nx (with x — A or B, the Fourier transform being 
represented by the ~ notation): 



Wx 



~1 2 

Px -Nx 



(7) 



As shown in Fig.|7| (6* Cen data), we estimate the PSD 
of the noise directly from the PSD of the signal by as- 
suming that the photon shot noise and detector noise are 
constant with respect to frequency (white noise). The con- 
tribution from the photometric fluctuations is visible at 
low frequencies. Considering the high frequency part of 
the spectrum, we extrapolate the level of the PSD back- 
ground to the lower frequencies. Due to the high SNR 
of the averaged PSD, the estimation of the background 
is precise enough to reconstruct Wa and Wb- The first 
frequency for which the signal becomes smaller than the 
noise marks the practical limit of the Wiener filter, and 
the higher frequency values are set to zero. This method 
is more efficient than estimating the noise PSD from the 
off-source batch, as it directly takes into account the pres- 
ence of photon shot noise, that also has to be removed. 
The resulting Wiener filters are presented in Fig.|S| The 
filtering of the photometric channels by these filters gives 
a clean P signal, as shown in Fig.lHl 



Fig. 8. Wiener filters computed from 9 Cen photometric 
signals (Fig.O. Wa is the upper Hne, and Wb the lower 
line. 
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Fig. 9. Photometric normalization signal P (thick line), 
with the Wiener filtered Pa (dashed line) and Pb channels 
(thin line). P is the geometric mean of the two filtered 
photometric channels. 



5.3. Alternative normalization methods 

If the SNR of the photometric channels Pa and Pb reaches 
too low values over the scan length, we choose to normalize 
the intcrferograms simply by averaging P over the fringe 
length, instead of using the Wiener filtered signal. This 
allows us to significantly reduce the amplification of the 
noise due to the normalization division. The limit between 
the two regimes is usually set to 5 times the readout noise. 
For intcrferograms that present very low photometric sig- 
nal over the fringe packet itself, we discard the scan as a 
significant bias can be expected on the modulated power. 
Both averaging and Wiener filtering are almost equiva- 
lent on the final calibrated intcrferograms, with a slight 
advantage to the Wiener filtering when the photometric 
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fluctuations are important (as in the UT observations for 
example) . 

5.4. Interferogram subtraction 

After their calibration, we subtract the two interferograms 
7i cai and I2 cai , in order to cancel the residual photomet- 
ric fluctuations due to the uncertainty in the estimation of 
the K coefficients. This subtraction has proven to signifi- 
cantly enhance the immunity of the interferograms to the 
contamination coming from the photometric fluctuation 
background. Fig. 1101 shows the calibrated and normalized 
interferometric signal hcai and l2cah together with / the 
result of the subtraction of these two signals defined as: 



2.5 ^ 



h. 



L 



2 cai 



(8) 



The combined signal / is used for the integration of the 
fringe power to derive the visibility (Sect.d). The advan- 
tage of using / instead of using separately Ii cai and I2 cai 
for the integration of the fringe power is that all the cor- 
related noise between the two signals is eliminated by the 
subtraction, while the fringes, perfectly opposed in phase, 
are amplified. This allows us to eliminate the residual pho- 
tometric fluctuations as well as part of the noise intro- 
duced during the photometric calibration. 

6. Interferogram selection 

6.1. Piston effect 

The photometric calibration of the interferograms com- 
pensates for the incidence of wavefront corrugation across 
each subpupil of the interferometer, however it does not 
help remove the random phase walk (differential piston) 
between the two subapcrtures. 

The differential piston, considered as a time-dependent 
OPD error x{t), can be locally expressed by a polynomial 
development around a reference time to (corresponding, 
for example, to the middle of the acquisition sequence): 



x{t) ^xo+x {t~to)+x {t- to)^ + 



(9) 



The effect of the OPD perturbation on the interferogram, 
and its consequence on the coherence factor measurement, 
depends on the order: 

— Zeroth order: the constant term xq can be seen as 
a global offset of the fringe packet. It is detected and 
corrected by the QL algorithm which centers the fringe 
packet in the middle of the interferogram. 

— First order: The first order of the piston x changes 
the fringe velocity and induces a simple frequency shift 
in the PSD. It modifies the fringe peak position, but 
acts only as a homothetic compression or expansion of 
the fringe packet along the OPD direction. The first 
order piston has no immediate effect on the fringe vis- 
ibility. However, if the shifting speed x of the fringe 
packet is too high, it can result in an undersampling 
of the fringes that will affect the visibility. 



1.5 
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0.0 — Afll 






20 40 60 80 100 

Optical path length (microns) 

Fig. 10. From top to bottom: /i normalized interferogram 
(black), I2 normalized interferogram (grey), and the result 
/ of the subtraction of these two signals. For clarity, the Ii 
and I2 signals are shifted vertically by -t-2 and -1-1, respec- 
tively. The correlated noise has disappeared in the com- 
bined signal and the fringe packet appears more clearly. 



— Second and higher orders: Any term of order two 
{acceleration) and beyond breaks the linear relation- 
ship in the scan between time and OPD, and conse- 
quently the Fourier relationship which is at the ba- 
sis of the visibility calculation, distorting the shape of 
the fringe peak. This introduces a non-linear, seeing 
induced multiplicative noise on the visibility measure- 
ments, which is the dominant noise source for strong 
signals (bright, unresolved objects). 

Detailed studies of the properties of atmospheric piston 
can be found for example in Linfield et al. H2(JU1|I and Di 
Folco et al. H2U(J2|I . When a dedicated fringe tracking in- 
strument becomes available on the VLTI (e.g. FINITO, 
Gai et al. 12003(1 . most of the piston will be actively re- 
moved by a servo loop. It remains to be checked how the 
residuals will still limit the final visibility precision. 
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6.2. Quality control of the interferograms 

The goal of the interferometric data processing is to ex- 
tract the squared visibility V^ of the fringes. The inter- 
mediate step to this end is to measure the squared coher- 
ence factor /i^ of the stellar light. This instrument depen- 
dent quantity characterizes the fraction of coherent light 
present in the total flux of the target. It is calibrated us- 
ing observations of a known star, as described in Sect.^ 
To avoid any bias on /i^, we have to reject the interfer- 
ograms that do not contain fringes (false detections), or 
whose fringes are severely corrupted by the atmospheric 
turbulence (photometrically or by the piston effect). The 
selection procedure is in practice similar to a shape recog- 
nition process. 

For this purpose, we measure in the wavelet power 
spectral density (WPSD) the properties of the fringe peak 
both in the time and frequency domains, and we subse- 
quently compare them with the expected properties of a 
reference interferogram of visibility unity, derived from the 
spectral transmission of the instrument. In this paper, we 
will refer interchangeably to the "time domain" or "optical 
path difference (OPD) domain" for the WPSD, as they are 
linearly related through the scanning speed of the VINCI 
piezo mirror that is used to modulate the OPD. The fringe 
peak is first localized in frequency by the maximum of the 
WPSD, and then the full width at half maximum is com- 
puted along the two directions: time and frequency. As the 
fringe packet has been recentered before the calibration, 
its position in the time domain is zero. Three parameters 
are then checked for quality: 

— peak width in the time domain (typically ±50% around 
the theoretical value is acceptable), 

— peak position in the frequency domain (±30%), 

— peak width in the frequency domain (±40%). 

In principle, the variation of the fringe contrast over the 
spectral band should also be taken into account to create 
the theoretical reference interferogram. But in practice, as 
long as the visibility of the fringes does not cancel out for 
a wavelength located inside the spectral band of the obser- 
vations, the shape of the interferogram remains very close 
to theoretical fringes of visibility unity. However, when the 
fringe visibility goes down to zero for a wavelength per- 
taining to the observation band, the fringe packet appears 
split in two parts in the time domain. Such a deviation 
from the "single packet" case can cause misidentifications, 
and eventually induce a bias on the derived /i^ value, as 
some valid interferograms would be rejected. In this case, 
one should use a dedicated broadband model to predict 
the reference interferogram, taking into account the ex- 
pected angular size of the target. The selection criteria 
should then be adapted to match this reference (increased 
peak width, presence of two maxima in the WPSD,...). 
An alternative is to directly adapt the basis functions of 
the wavelet transform so that they match the predicted 
interferograms. We have not implemented these methods 
in our current processing algorithm, whose validity is thus 
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Fig. 11. Wavelet power spectral densities of processed 
interferograms {a represents the wave number, i.e. the 
inverse of the wavelength). The upper figure shows the 
WPSD of a good quality interferogram: the energy is well 
confied in the fringe power peak. The bottom WPSD is 
affected by strong atmospheric piston: a significant part of 
the fringe power is spread outside of the theoretical fringe 
peak, both in time and frequency domains. The isocon- 
tours delimit 77 % and 76 % of the total modulated power, 
respectively. 



limited to the cases when the visibility is above zero for 
all wavelengths in the observation band. If this condition 
is not realized, then the interferogram selection should 
be disabled, or at least the selection criteria should be 
made significantly less stringent, in particular regarding 
the fringe peak width. 

Fi ellll shows two examples of interferogram WPSD, 
one of them being affected by atmospheric piston. The 
difference in terms of fringe peak shape is clearly notice- 
able, and leads to the rejection of the corrupted interfer- 
ogram (bottom figure). This selection process has shown 
a very low false detection rate, and rejects efficiently the 
interferograms that are affected by a strong piston effect. 
However, limited piston of order two (and above) is not 
identified efficiently. The problem here is that the relevant 
properties for the estimation of the second order piston are 
currently difficult to measure with a sufficient SNR from 
the data, as they are masked by the order 1 piston. We ex- 
pect that the introduction of the FINITO fringe tracker in 
2004 will allow us to derive an efficient metric to reject the 
interferograms affected by a high order piston effect. After 
the fringe power integration (described in Sect.[7I), we fil- 
ter out the scans which //^ deviates by more than 3 a from 
the median of the full batch of interferograms (usually 500 
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scans). This step prevents the presence of very strong out- 
hers, which can appear due to the division introduced by 
the normaUzation of the interferograms (introduction of 
Cauchy statistics). 



6.3. Immunity to selection biases 



An essential aspect of the parameters used for the quality 
control of the fringe peak properties is that they are largely 
independent of the visibility of the fringes, and therefore 
do not create selection biases. In particular, the integral 
of the fringe peak (directly linked to the visibility) or its 
height are never considered in the selection. The parame- 
ters chosen in Sect. l6.2l clearlv depend on the photometric 
SNR, but are independent of the visibility of the fringes, 
thanks to the calibration procedure described in Sect.|Sl 

The upper part of Fig. II 21 demonstrates this indepen- 
dence in the difficult case of the batch of interferograms 
obtained on a Cen A. Despite the very low visibility of the 
fringes, no systematic deviation is visible for low photo- 
metric SNR values, as the dispersion is symmetric around 
the mean value. The same plot for 9 Cen f Fig. 1121 bottom) 
does not show any deviation either. A further discussion 
of the properties of the histogram of these measurements 
can be found in Sect. 18.31 This means that the quality 
control described in this paragraph is not linked to the 
observable, and thus does not introduce a selection bias. 
Its effect in the case of the 9 Cen and a Cen A observations 
is discussed in Sect. 18.31 

A critical case is when the visibility is extremely low. 
In this situation, the fringe peak will tend to blend in 
with the noise, which tends to make it appear broader 
and slightly displaced. Therefore, low-visibility data are 
more likely to be rejected than high-visibility data. This 
can introduce a bias towards higher /i^ for low-visibility 
observations: a scan with a -|-1 ct deviation is accepted, but 
a scan with a -1 cr deviation is more likely to be rejected 
as it fails the selection criteria. However, in this situation, 
the risk is high to fail to reject the spurious spikes that 
are created in the calibrated interferograms due to the 
division by the P signal (Sect.EJ. Without the selection 
procedure, the modulated power of these calibration arte- 
facts will be integrated in the final /i^ value. As this power 
is essentially random, but always positive, these misiden- 
tifications would then result in a strong positive bias on 
the final ^^ value. For this reason, and in spite of the 
potential rejection of a small part of the valid interfero- 
grams, the application of the selection procedure results 
in a more reliable estimate of /i^ , even for the very low vis- 
ibility fringes. In any case, the careful examination of the 
statistical properties of the mu^ histogram (see Sect. |H1), 
and in particular of its skewness, allows us to detect a 
possible selection bias. 
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Fig. 12. Squared coherence factor /x^ as a function of the 
photometric signal on the fringe packet, for the sample 
batches obtained on a Cen A (top) and 6* Cen (bottom). 
The dashed lines show the mean value of the series derived 
using the bootstrap technique. 



7. Visibility computation 

7.1. Wavelet power spectral density 

Once the interferogram has been calibrated and normal- 
ized, the squared coherence factor /i^ is measurable as 
the average modulated energy of the interference fringes 
over the batch. It is computed by integrating the power 
peak of the interferograms in the average WPSD (see also 
Appendix El and Sect. l7.2.1|) . The WPSD is a two dimen- 
sional matrix, examples of which are shown in Fig llll For 
all wavelet transforms, we use the Morlet wavelet, which 
is defined as a plane wave multiplied by a Gaussian en- 
velope. It closely matches the shape of the interferomet- 
ric fringe packet. When computing a classical PSD, the 
interferogram is projected on a base of sine and cosine 
functions, which are not localized temporally. This means 
that the information of the position of the fringe packet 
is not used, and that the noise of the complete interfer- 
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ogram contributes to the measured power. On the other 
hand, the wavelet transformation projects the interfero- 
gram on a base of wavelets that are localized both in time 
and frequency, making full use of the localized nature of 
the modulated energy. As discussed in Appendix ^ the 
modulated energy of the signal is conserved by the wavelet 
transform in the same way as through the classical Fourier 
transform. 

7.2. Estimation of the squared coherence factor jj.'^ 

7.2.1. Fringe power integration 

The average power spectral density of the a Cen A sample 
batch, computed using the wavelet transform, is shown 
in Fig.El To obtain this ID spectrum from the original 
2D WPSD matrices, we first project the WPSD matrix of 
each interferogram on the frequency axis, by integrating 
it over the fringe packet length (time axis). From this we 
obtain a series of one-dimensional vector PSDs, similar to 
the Fourier PSD but with a reduced noise. Before the av- 
eraging, we recenter each fringe peak using the frequency 
position information derived from the selection of the in- 
terferograms (Sect-EJ. This step allows us to confine more 
tightly the energy of the peak, which is displaced by the 
first order piston effect. This reduces the influence of pis- 
ton on the final /z^ value. The co-added ID spectrum is 
the signal used for the final power integration to estimate 
the //^ of the star. 

The integration of the fringe power is typically done 
over 100 pixels in the time domain (20 fringes), and from 
2 000 to SOOOcm^^ in the frequency domain (see also 
Sect. EH) • 

7.2.2. Removal of the WPSD background 

The power in the fringe peak is contaminated by three 
additive components: 

— the photon shot noise, 

— the detector noise, 

— the residual photometric fluctuations. 

To estimate the modulated power of the fringes, it is es- 
sential to precisely remove these contributions from the 
PSD of the interferograms. 

Perrin (2003:1 has developed an analytical treatment of 
the photon shot noise based on its particular properties 
(Poisson statistics). The photon shot noise is perfectly 
white, as it is created from a purely random process. 
However, due to the calibration and normalization pro- 
cess of the interferograms, its translation onto the final 
interferometric signal / could in theory deviate from this 
property and show a dependence with frequency. Such an 
effect has not been observed in practice on the VINCI 
data, and the uniform subtraction of the photon noise 
background from the PSD of the / signal has proven 
to be very efficient. A good example of the "whiteness" 
of the photon shot noise of the processed fringes can be 
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Fig. 13. Average WPSD of the a Cen A observation. The 
original two-dimensional matrix has been integrated over 
the fringe packet length in the time domain. The resulting 
projection on the wave number axis allows one to visualize 
clearly the noise contribution. The subtracted noise model 
is shown by the dashed line. The final WPSD (thick line) 
shows no bias in spite of the large brightness of the star 
and very low squared coherence factor of the fringes {/j,^ ~ 
0.3%). 



found in Wittkowski et al. H2()()4|l . where a very bright star 
{rriK = —0.6) was observed with the two 8 m telescopes 
UTl and UT3 {B = 102.5m). In spite of the extremely 
large flux on the VINCI detector (100 m^ collecting optics) 
and the very low visibility of the fringes {V"^ ~ 10~^), the 
resulting PSD background is white, therefore validating 
our photon shot noise removal method under the most 
demanding conditions. 

In order to fully justify our background removal proce- 
dure, we still have to verify the "whiteness" of the detector 
noise, whose statistics and frequency structure depends on 
the type of detector and readout electronics used. The in- 
frared camera of VINCI (LISA) is based on a HAWAII 
array, which is read using an IRACE controller (Meyer et 
al. i2Q00) . As only four pixels of the 1024x1024 array are 
actually used, an engineering grade detector was chosen 
for the instrument. It presents a large quantity of dead 
and hot pixels, and therefore it was necessary to thor- 
oughly check its noise characteristics. This was achieved 
during extensive laboratory testing, and is also verified 
automatically for each observation. It appeared that the 
LISA detector noise is perfectly white, without any signif- 
icant electronic interference signature. 

This satisfactory behavior of the detector and photon 
shot noises allows us to remove them simultaneously by 
subtracting to the ^^ value of each interferogram an aver- 
age of its WPSD at high frequency, measured outside of 
the domain of frequency of the interferometric fringes. To 
correct for potential residuals pf the photometric calibra- 
tion, we fit a linear model of the residual background to 
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the average WPSD of the interferograms in the batch. In 
this procedure, we aUow for a Hmited slope of the back- 
ground model, in order to correct a possible residual power 
from photometry. Thanks to the averaging of a large num- 
ber of scans, the noise on the average WPSD low, and 
the fitting procedure is very precise. Most of the time, 
and even for the most important fluctuation cases (Unit 
Telescopes in multispeckle mode) , the contribution of the 
residual photometric noise is totally negligible on the com- 
bined interferogram / obtained from the subtraction of the 
calibrated signals Iicai and l2cai (Eq. |SJ). 

An illustration of the background quality is presented 
in Fig. El The WPSD background noise appears perfectly 
white, even at the very enlarged scale used to visualize 
the very small fringe peak of aCenA (/i^ « 0.3%). In 
particular, no "color" or electronic interference ("pickup") 
are present. The model of the removed noise is also plotted 
on the two figures, showing that the subtraction is very 
clean. 



7.3. Estimation of the statistical error 

To compute the statistical error on the /i^ estimation, we 
integrate separately the fringe power in each WPSD of the 
batch, correct the detector and photon shot noise biases 
individually, and use a weighted bootstrapping technique 
on this set of measurements. Our sample is made of N 
pairs {isf,Wi) where /zf is the squared coherence factor 
obtained by integrating the WPSD of the scan of rank i 
in the series and Wi is its associated weight. It is defined 
as the average level of the photometric signal P over the 
fringe packet length (20 fringes in the K band) multiplied 
by the inbalance between the two photometric channels 
Pa and Pb'. 



mmiPA.t,PB.ij I fT, — 5— 

' Fringes 



max{PA,i,PB,i) 



(10) 



It characterizes well the clarity of the total photometric 
signal that contributes to the formation of the fringes. 
The final dispersion of the ji^ values is reduced by this 
weighting. The detailed description of the weighted boot- 
strapping method used for the computation of the error 
bars is given in Appendix IbI 

The bootstrapping technique has the important advan- 
tage of not making any assumption on the type of statisti- 
cal distribution that the data points follow. In particular, 
it is more reliable than the classical approach that assumes 
a Gaussian distribution of the measurements. Skewness 
and other deviations from a Gaussian distribution are au- 
tomatically included in the error bars, which can be asym- 
metric. 

The statistical dispersion of VINCI measurements 
shows two regimes: for bright stars the precision is lim- 
ited by the piston and photon shot noise, while for the 
fainter objects, the main contributor to the dispersion is 
the detector noise of the camera, and the precision de- 
grades rapidly. A discussion of the different types of noise 



intervening in the visibility measurements can be found 
for instance in Colavita 1)1999(1 and Perrin 1)2003(1 . The p? 
measurements discussed in this paper have a relative sta- 
tistical precision of ±3.00% for aCenA, and ±0.53% for 
0Cen. The lowest relative statistical dispersion cr(/z)/// 
reached up to now on the coherence factor with VINCI 
is in the 2% range. Under good conditions, this translates 
into a bootstrapped statistical error of less than ±0.1% on 
/i for 5 minutes of observations. 



8. Post-processing quality control 

After a batch of interferograms is processed, several qual- 
ity controls are performed in order to detect any prob- 
lem in the resulting visibility values and statistical error 
bars. This step is essential to ascertain the quality of the 
interferometric data, as it can vary depending on the at- 
mospheric conditions (e.g. seeing, coherence time) and on 
the general behavior of the instrument (e.g. injection of 
the stellar light in the optical fibers, beam combiner prop- 
erties, polarization mismatch of the two beams). 



8.1. Power peak integration boundaries 

A potentially damaging effect of the atmospheric piston 
on the visibility of the fringes is that it tends to move 
the position of the fringe peak, and to spread it over a 
wider frequency range. If the frequency boundaries for 
the integration of the fringe peak are set too tight, the 
result could be that part of the modulated power is not 
taken into account, creating a bias. These boundaries are 
automatically changed as a function of the ground base- 
line length to account for the increased piston strength on 
longer baselines. They are not modified as a function of 
the projected baseline, and are thus identical for scientific 
targets and calibrators. 

To check for the presence of such an effect, we mea- 
sure the fringe peak shape in the WPSD. More precisely, 
we estimate its central wave number, full width at half 
maximum, as well as the limit wave numbers for which 
the background level is reached. Using these extended lim- 
its, we integrate the fringe power and compare this value 
to the one obtained with the user-specified wave number 
limits. If a discrepancy is found at a significant level, the 
batch is considered dubious and can be rejected after fur- 
ther examination. 



8.2. Histogram properties 

As the noise sources acting on the /i^ values have normal 
statistics, it is expected that the distribution of the /x^ 
values over the batch is also normal. Although the boot- 
strapping procedure used to compute the /i^ error bars is 
not sensitive to the type of distribution, a large skewness 
or kurtosis would betray a problem in the cahbration of 
the interferograms that could eventually bias the final ji^ 
value. The relevant parameters for this verification are the 
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Table 4. Reasons for the rejection of 9 Cen and a Cen A 
interferograms during the processing. The lower part of 
this table corresponds to the selection criteria related to 
the atmospheric piston effect. 

Reason 

Photometry too low 
Large OPD jump 
Fringes at edge 
Fringe packet width 
Fringe peak position 
Fringe peak width 
Outliers (3 a) 

Total number of rejected scans 114/500 223/500 



skewness coefficient s (third moment of the distribution) 
defined as: 



Cen 


aCen A 


77 


24 


13 


47 


6 


27 


1 


47 


3 


40 


5 


33 


9 


5 






l)a3 



and the kurtosis coefficient k (fourth moment): 



N 



2^4 



^ (iV - 1)(t4 



(11) 



(12) 



where ^f are the squared coherence factor values, a the 
standard deviation, /i^ the unweighted average, and N the 
number of scans. The skewness characterizes the presence 
of a "tail" on the histogram. A large value of s is therefore 
a symptom of a potential bias problem in the distribution, 
as a significant number of values are either too large or too 
small compared to the average value of the sample. A pos- 
itive kurtosis coefficient means a distribution more peaked 
than the normal one. However, it should be stressed that 
the kurtosis is not a very robust parameter to assess if the 
sample is drawn from a normal distribution. R requires 
a large number of sample values to be relevant, and it 
is very sensitive to the presence of outliers. Therefore, it 
should only be used in conjunction with other statistical 
tests. A range of ±0.5 can be considered acceptable for 
k. When random samples are drawn from a normal pop- 
ulation, the resulting skewness coefhcients will fall into 
the range ±0.18, with a probability of 90%. We therefore 
choose this value as an acceptable range. 

8.3. Application to 9 Cen and a Cen A 

Table 21 gives the reasons for the rejection of the interfer- 
ograms of the 9 Cen and a Cen A batches. In the case of 
a Cen A, a larger number of interferograms are rejected 
due to the very low visibility of the fringes. 

The measured statistical properties of the processed 
interferograms of 9 Cen and a Cen A are given in Table [S] 
The values in brackets were obtained by disabling the pis- 
ton selection of the interferograms (based on the fringe 
packet width, and on the position and width of the fringe 
peak in the power spectrum of the interferogram) . The 



Table 5. Statistical properties of the 6'Cen and a Cen A 
sample batches. The values obtained when disabling the 
selection of the interferograms based on their piston prop- 
erties are given in brackets for comparison. 



eCen 



a Cen A 



Reduced scans 
Average fj,^ (%) 
Stat, error (la) 
Rel. error a/fi^ 
Skewness s 
Kurtosis k 



386 (395) 

8.995 (8.999) 

0.048 (0.050) 

0.53% (0.56%) 

-H0.023 (+0.013) 

+0.164 (+0.062) 



277 (391) 
0.3180 (0.2949) 
0.0095 (0.0142) 
3.00% (4.81%) 
-0.007 (+0.160) 
+0.044 (-0.042) 
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Fig. 14. Histogram of the fj,^ values obtained on 6 Cen. 
No significant skewness is present. 



comparison of the selected vs. non-selected versions of the 
data processing shows that the piston selection has a pos- 
itive effect on the dispersion of the measurements. For 
9 Cen, the difference is minimal between the two kinds 
of processing. In particular, the total number of pro- 
cessed interferograms is almost identical for the two cases. 
However, for a Cen A, the difference is clearly noticeable, 
as the final error bars are 60% larger when the selection 
is disabled, in spite of a total number of processed scans 
approximately 40% larger. The skewness of the histogram 
is also much larger in this case (by a factor of 20). This 
clearly shows the advantage of the fringe selection pro- 
cedure, in particular for the rejection of the calibration 
artefacts (false detections) in the very low visibility case 
(see also Sect. 16.3(1 . 

In the case of 6'Cen (Fig.[l3|) and a Cen A (Fie-. IT^ . 
no skewness is present. For 9 Cen, a small positive kurto- 
sis fc « 0.16 is detected, meaning that the distribution is 
slightly too peaked {leptokurtic, as opposed to a platykur- 
tic distribution that is too flat). However, it is easily inside 
the acceptable range (±0.5), and this property is taken 
into account in the bootstrapped error bars. 

9. Visibility calibration 

9.1. Principle 

The data reduction software of VINCI yields accurate es- 
timates of the squared modulus of the coherence factor 
fi'^, which is linked to the squared object visibility V"^ by 
the relation: 



y2^^ 



rp2 



(13) 
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Fig. 15. Histogram of the fi^ values obtained on a Cen A. 
A moderate skewness is present towards higher ^^ values. 



where T^ is the response of the system to a point source, 
also called transfer function (hereafter TF), interferomet- 
ric efficiency, or system visibility. It is measured by brack- 
eting the science target with observations of calibrator 
stars whose V'^ is supposed to be known a priori. The 
accuracy of our knowledge of the calibrator angular diam- 
eter, and the precision with which we estimate its /i^ are 
therefore decisive for the final quality of the scientific tar- 
get observation. Typically, the scientific targets are brack- 
eted by calibrator observations, so as to be able to verify 
the stability of the TF. In this respect, the VLTI has often 
proved to be stable at a scale of a few percent over several 
nights. Nevertheless, and to guarantee the quality of the 
VINCI data, calibrators are regularly observed during the 
night, before and after each scientific target. 

9.2. Interferometric transfer function estimation 

9.2.1. Transfer function model 

By nature, the interferometric TF is affected by a large 
number of parameters: atmospheric conditions (seeing, co- 
herence time) , polarization (incidence of the stellar beams 
on the siderostat mirrors, spectrum of the target, etc.). 
These effects combine to make T^ a stochastic variable, 
that can evolve over a wide range of timescales. In or- 
der to estimate its value and uncertainty on a particular 
date at which it was not directly measured (e.g. during 
the observation of a scientific target), it is necessary to 
use a model of its evolution. Such a model relies necessar- 
ily on an hypothesis, for instance that the value of T^ is 
constant between two (or more) calibrator observations, 
that it varies linearily, quadratically, or any higher order 
model. Let us now evaluate the most suitable type of TF 
model for the observations with VINCI. 

As a practical example. Fig. Ilfilshows the evolution of 
T^ over one night of observations, with a typical sampling 
rate of one measurement every 15 minutes. This series of 
27 observations was obtained during the night of 29 May 
2003 on the EO-GO baseline (16m ground length). A num- 
ber of different stars with known angular diameters were 
observed, covering spectral types in the G-K range. During 
these observations, spread over 8 hours, the seeing evolved 
from 1.0 to 2.0 arcsec, the altitude of the observed objects 
was distributed almost uniformly between 25 and 80 de- 
grees, and the azimuth values covered 15 to 90 degrees (N 





















- 




















: 


^ 42.0% - 






n T 


11 


-I. 




I 


- 




■i 


r 


% 


. 






■ 




If 












39.0% - 





















3.0 4.0 5.0 

Relative time (h) 



Fig. 16. Evolution of the transfer function T^ during one 
night (2002-05-29) on the EO-GO baseline of the VLTI 
(16m in ground length). Each symbol corresponds to a 
different star. 



= 0, E = 90). Due to this broad range of conditions, this 
series represents a worst case in terms of TF stability. As 
a reminder, under normal conditions, a calibrator is se- 
lected as close as possible to the scientific target, in time, 
position and spectral type. 

Over the whole night, the overall stability is satisfac- 
tory, with a dispersion of CTtot — 0.64% around the average 
value of T^ = 41.75%. In order to estimate the external 
dispersion CTcxt of the transfer function over the night (due 
to the atmosphere and instrumental drifts), we can sub- 
tract the average of the intrinsic variances af of the T^ 
values crfjjj from the total variance a^^^: 

c^oxt = f^tot - o-fnf (14) 

The average precision of each individual T^ measurement 
in our sample night is CTint = 0.21%. This gives an ex- 
ternal dispersion of (Text — 0.60%. In this particular case, 
the external dispersion is thus dominant over the internal 
measurement errors, by a factor of almost three. 

From this example, we can conclude that the rate of 
one measurement every 15 minutes is insufficient to sam- 
ple the fluctuations of the TF. Due to this, we do not gain 
in precision by interpolating the TF values using a high 
order model (quadratic, splines,...). In the current state 
of the VLTI (siderostat observations), the most adequate 
model for the estimation of the TF is thus a constant value 
between the observations of the calibrators. The 1.8 m 
Auxiliary Telescopes will soon allow us to sample the TF 
with a much higher rate, of the order of 1 minute, and 
higher order models of the TF variations could become 
necessary. As we are dominated by the external dispersion 
(Text, the uncertainty on the TF has to be estimated from 
the dispersion of the individual T^ measurements obtained 
before and after the scientific target, without averaging of 
their associated error bars. 
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9.2.2. Refinement of the transfer function model 

Under good and stable conditions, the random disper- 
sion of T^ introduced by the atmosphere can be very low 
between two consecutive observations of a calibrator. In 
this case, we want to evaluate the true uncertainty on the 
model T^ by comparing the hypothesis of stability to the 
calibrator observations, and subsequently refine the hy- 
pothesis used to estimate the error bar on T^. 

The observational strategy chosen with VINCI is to 
record several series of interferograms consecutively for 
each calibrator observation (typically three), over a pe- 
riod of about 15 minutes. To decide if the atmospheric 
and instrumental conditions are stable over this period, 
we compute the following x^ expression: 



Xtot 



E 



(Ti-T^) 



cr. 



(15) 



stat,z 



where Tf are the consecutive estimates of the TF obtained 
on the calibrator, a'^^.^^ ^ the statistical error of each mea- 
surement, and T^ the weighted average of the Tf values 
(using the inverse of the statistical variance as weights). If 
the resulting y^ is small (less than 3), then the hypothesis 
that the TF is constant is probably true: the Tf values 
can be averaged and the global statistical error bar re- 
duced accordingly. If not, then this hypothesis cannot be 
made, and a realistic approach is to consider as the true 
measurement error of the average TF the standard devia- 
tion of the Tf sample. 

When several series of interferograms are obtained on 
the same calibrator and the conditions described above 
are verified, the resulting estimates of the TF can be av- 
eraged in order to reduce the attached statistical error bar. 
However, the systematic error introduced by the a priori 
uncertainty on the angular size of one calibrator cannot 
be reduced by repeatedly observing this star, but only by 
combining the TF measurements obtained on independent 
objects. 



9.3. Application to the sample observation of a Cen A 

In the case of the observations described in this paper, 
9 Cen was observed one hour before a Cen A. Assuming a 
UD angular diameter of 6'ud ~ 5.305±0.020mas (Kervella 
et al. |2£)03b), and taking into account the spectrum of 
the source and the bandwidth averaging effect (also called 
bandwidth smearing, see e.g. Davis et al. l2()()()|l . we expect 
a squared visibility of V^^oo ~ 0.1796 and a 1 cr systematic 
uncertainty agyst = ±0.0027 for the 65.929 m projected 
baseline (weighted average over the interferogram series) . 
As we observed an instrumental coherence factor of fi^ = 
0.08995 ± 0.00048, the transfer function T^ is estimated 
to be: 



^2 _ 0-08995 [±0.00048]stat 
0.1796 [±0.0027]syst 



Using the relations detailed in Appendix El to compute 
the error bars of the ratio fM^ /T'^, we obtain 



t2 = 0.5010 [±0.0027],tat [±0.0077],y,f 



(17) 



As we consider here only one calibrator observation, we 
cannot estimate the external dispersion of T^, and we 
consider only the internal statistical and systematic error 
bars. As a remark, this T^ value is not identical to the one 
computed in Sect. l9.2~n as it was obtained more than one 
month later. The VINCI coupler is known to be sensitive 
to long term temperature variations (over a timescale of 
weeks) , an effect that can explain the observed difference. 
The squared visibility value of a Cen A is then: 



l^aCcn = 0.00635 [±0.00019]stat [±0.00010]syst. 



(18) 



(16) 



The uncertainty on this value is dominated by the statis- 
tical error, despite the importance of the systematic error 
on the value of T^. The average baseline of this measure- 
ment is 61.470 m, and we can now deduce the uniform 
disk model angular diameter of a Cen A, 6'ud = 8.305 ± 
0.024 mas, which is very close to the published value of 
6'ud = 8.314 ± 0.016 mas from Kervella et al. H2003b|l . 
This computation takes into account the wavelength av- 
eraging effect due to the broadband K filter of VINCI as 
described by Kervella et al. I)2003b|l . 

10. Conclusion 

We have described the data reduction methods that are 
used on VINCI, the VLTI Commissioning Instrument. In 
particular, we detailed the photometric calibration of the 
interferometric signals, followed by the normalization of 
the fringes, and the subtraction of the two calibrated 
interferograms. Due to the efficient spatial filtering pro- 
vided by the single-mode optical fibers, this procedure 
provides a clean calibration of the fringes, and allows us 
to derive the squared coherence factor /i^ with high accu- 
racy. Combined with observations from a calibrator star, 
it yields the squared visibility V'^ . This value can be inter- 
preted physically through the use of a dedicated model of 
the observed object. Applying the data reduction methods 
described in this paper to sample data from a Cen A yields 
a realistic value of its uniform disk angular diameter. Our 
procedures can easily be adapted to other single mode in- 
terferometric instruments. In particular, they can be gen- 
eralized to spectrally dispersed fringes and to a multiple 
beam recombiner using the integrated optics technology 
(Kern et al. 2003 ) . Such a device could allow the simulta- 
neous recombination of the beams from the four 8 m Unit 
Telescopes and four Auxiliary Telescope of the VLTI in a 
compact instrument. 

Appendix A: The wavelet transform 

The wavelet transform belongs to the class of time- 
frequency transforms which are powerful tools to study 
non-stationary processes such as turbulent flows in fluid 
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mechanics. Wide band coaxial interferogranis recorded 
through a turbulent atmosphere can be strongly distorted 
due to the differential piston effect and fast photomet- 
ric fluctuations. In this context, the wavelet transform is 
an efficient tool to study and analyse the interferograms 
recorded from the ground. 

The continuous wavelet transform (hereafter CWT) is 
defined by: 



W{s,t) 



1 



f{t)r 



dt 



(A.l) 



where f{t) is the signal defined as a function of time, </> (i) 
the chosen wavelet function, "0* its complex conjugate, s 
the scale, and r the translation. 

For the present application of the CWT to interfer- 
ometry, we have chosen to use the Morlet wavelet, that is 
defined as a Gaussian envelope multiplied by a plane wave 
(Goupillaud et al. 1984 , Farge 1992) : 



V' (??) = exp {i2'KVaT]) exp (—77^/2) 



(A.2) 



where ry is the non dimensional time parameter and ui^ 
is the non dimensional frequency. Initially used for the 
analysis of seismic signals, the Morlet wavelet is a good 
approximation of the fringe pattern produced by VINCI. 
The data processing methods presented here make exten- 
sive use of this particular wavelet for the recognition and 
localization of the fringes (Sect. ^3 ^^^d subsequently for 
the integration of the modulated power of the interference 
fringes (Sect.O. Fig lA.ll shows the shape of the imaginary 
part of the Morlet wavelet assuming typical parameters for 
the processing of data from the MONA beam combiner {K 
band). 

If we now express the CWT in the Fourier domain 
(Eq. IA.3|I , it appears clearly that the CWT is a filtered 
version of the signal for different sets of filters: 



W{s,t) 



+ 00 



/ [v) r M e' 



dv 



(A.3) 



Since the CWT is simply a convolution between the sig- 
nal / {t) and expanded/contracted versions of the wavelet 
function (Eq. lA.ip , the Morlet wavelet is very efficient to 
analyse wide-band coxial interferograms. The CWT of an 
interferogram using the Morlet wavelet is a complex quan- 
tity and its maximum energy is found for the wavelet that 
is most similar to the recorded interferogram. 

The CWT using the Morlet wavelet is not orthogo- 
nal but since it relies on a set of filtered versions of the 
signal with strong redundancy, the original signal can eas- 
ily be reconstructed fFarge [T992I Perrier [l995|l . The en- 
ergy properties of Wavelets are similar to the ones of the 
Fourier analysis, with the equivalent of the Parseval the- 
orem fPerrier fl995|l . We have therefore the equivalence of 
the two following expressions of the energy E of the signal: 



E 



2G 



+ 00 /' + 00 



i) Jo 



\W{s,Tf%dT 



(A.4) 
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Fig. A.l. VINCI interferometric fringes (upper curve, 
from a processed interferogram of 6 Cen) and the Morlet 
wavelet function imaginary part (bottom curve). 



E 



+ 00 



/M 



dv 



with the coefficient C,i, defined as 



C^ 



+00 



i^(5^)r 



^ 



dv 



(A.5) 



(A.6) 



As a consequence, we are able to recover the modulated 
energy of the original signal (the interferometric fringes) 
by integrating its wavelet power spectrum over the time 
and frequency regions where the interferogram is present. 
Compared to the classical Fourier analysis, such an 
approach allows to minimize the biases due to both the 
white and colored (frequency dependent) noises. Thanks 
to its localization both in time and frequency, the Morlet 
wavelet is better suited to the study of interferometric 
fringe packets than the classical Fourier base functions 
(sine and cosine), as the noise present outside of the 
fringe packet in the scan is excluded from the integrated 
power. The interested reader will find a more detailed 
treatment of the wavelet transform in Daubcchies H1992|l . 
Farge jfTM^ .Perrier ifBMI and Mallat (1999) . 

Appendix B: Computation of weighted 
bootstrapped error bars 



Originally developed by Efron (1979^, the bootstrap anal- 
ysis, also called sampling with replacement, consists of 
constructing a hypothetical, large population derived from 
the original measurements and estimate the statistical 
properties from this population. This technique allows us 
to recover the original distribution characteristics without 
any assumption on the properties of the underlying pop- 
ulation (e.g. Gaussianity) . An introduction to bootstrap 
analysis can be found in Efron 1)1993(1 and Babu H1996(l . 

Our implementation of the bootstrapping technique 
draws, with repetition, a large number M of samples 
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containing N elements from the original set of measure- 
ments {^^,Wi), also A'' elements in length. ^| designates 
the squared coherence factor associated with the scan of 
rank i in the series, and Wi is its associated weight. The 
result of this drawing is an A^ x M matrix of {nl j,Wk.j) 
pairs (1 < fc < iV; 1 < i < M). The fact that the same 
element of the original data set can be repeated several 
times in the drawing is essential, as it allows us to create 
independent samples. Typically, several thousand samples 
are obtained from the original data, which contains a few 
hundred ^^ values. The weighted average values /z^ are 
computed for each of the N drawings: 



Ml 



1 



N 

Lfc=i wk,j fe=i 



(B.l) 



The resulting ensemble [/i^] (M elements) is sorted in as- 
cending order, and reindexed with the percentiles of the 
rank of each value in the set: 



2 
Mo/M ' 



2 



2 

' ^M/Af 



(B.2) 



The 16% lower and upper values are discarded, and the 
new extremes values of this vector give the limits of the 
68 % confidence interval: 



2 _ 2 . 2 _ 2 

Mmin — M 0.16 > Mmax — M 0.84 



(B.3) 



This is the probabilistic definition of 1 cr error bars, and 
we therefore obtain the <t^ and (t_ asymmetric error bars 
through: 



Ml, 



M^ 



M" 



2 
Mm in 



(B.4) 



where /i^ is the weighted average of the original sample 
[/if] using the weights [wi]. The same process can be ap- 
plied using 2.5-97.5% percentile limits to obtain the error 
bars equivalent to 2 a, and 0.5-99.5% for 3 a. 

Alternatively, one can derive the bootstrapped vari- 
ance (Tgg directly from the /i^ ensemble: 



1 *^ 

-is = ^ E (as' - 



(B.5) 



j=i 



The internal bias fees of the population is given by: 
^bs-mJ-A? (B.6) 

This bias is naturally taken into account in the computa- 
tion of the confidence interval limits as described above. 

Appendix C: Statistical and systematic errors of 
the ratio of //^ and T^ 

In the expression of T^ of Eq.^1 we have to separate the 
contributions from the systematic uncertainty on the cal- 
ibrator knowledge, and the statistical error of the instru- 
mental measurement of fj,^. While these two terms con- 
tribute to the global uncertainty on the squared visibility 
V'^, their nature is fundamentally different. While it is 



possible to reduce the statistical error by averaging sev- 
eral measurements, the systematic uncertainty originating 
in the calibrator diameter error bar will not be changed. 
This last term is therefore a fundamental limitation to 
the absolute precision of the visibility measurement. This 
limit can be reduced by using several calibrators, or by 
selecting very small stars as calibrators. We then benefit 
from the fact that the visibility function V^{B,9) for a 
stellar disk is nearly flat when the star is not significantly 
resolved, and the resulting systematic error on V^ remains 
small. 

Considering a symmetric error bar on the assumed an- 
gular diameter of the calibrator, the resulting error bar on 
the V^ estimate is in general not exactly symmetric, due 
to the non linearity of the visibility function. In practice, 
asymmetric error bars are easily manageable numerically. 
However, in order to simplify the notations in the present 
discussion, we make the assumption that this asymmetry 
is negligible. 

The estimation of the two kinds of error contributions 
relies on an approximation of the Cauchy statistical law 
characteristics. When dividing two normal statistical vari- 
ables X and y of respective means and standard devia- 
tions {x,a'^) and (y, ct^), the resulting ratio x/y follows a 
Cauchy statistics that has, strictly speaking, no defined 
mean value. It is therefore necessary to make an approx- 
imation for the case when ay <C y. In this case, a second 
order approximation of the mean 'z and variance cr^ of 
z = x/y is given by Browne H2002|l : 



(C.l) 




y 



4+S 



(C.2) 



It is not possible to obtain a meaningful average value of 
the ratio x/y if the standard deviation ay of the denom- 
inator y is not small compared to its average value. As 
a remark, the average value of x/y is in general different 
from x/y. 

The average value of the transfer function T^ and its 
associated statistical error bars are computed by replacing 
in formulas IC.ll and IC.2I the values of (x, u^) and {y,ay) 
by the following terms: 
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(C.3) 
(C.4) 



Similarly, the systematic error is computed using the re- 
placements: 



V 



A* JeCcn 
2 l2 



K 



thooJ Cen 







[<^thco\e Ceil 



(C.5) 
(C.6) 



Applying this computation to the numerical values found 
for 6 Cen, we find: 



T^ = 0.5009 [±0.0027]stat [±0.0077] 



syst 



(C.7) 
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We obtain the cahbrated squared visibihty of a Cen A: 
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